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Figure of Merit for Dark Energy Constraints from Current Observational Data 



OO 
O 

o 

(N 



m 
si 

6 



> 

in 
o 

oo 
O 



X 



Yun Wang 

Homer L. Dodge Department of Physics & Astronomy, Univ. of Oklahoma, 
440 W Brooks St., Norman, OK 73019; email: wang@nhn.ou.edu 

In order to make useful comparisons of different dark energy experiments, it is important to 
choose the appropriate figure of merit (FoM) for dark energy constraints. Here we show that for a 
set of dark energy parameters {/,}, it is most intuitive to define FoM=l/ ^/det Cov(/i, fa, fa, ■■■), 
where Cov(/i, fa, fa, ...) is the covariance matrix of {/»}. In order for this FoM to represent the 
dark energy constraints in an optimal manner, the dark energy parameters {fa} should have clear 
physical meaning, and be minimally correlated. We demonstrate two useful choices of {fa} using 182 
SNe la (from the HST/GOODS program, the first year Supernova Legacy Survey, and nearby SN la 
surveys), [R(z t ), l a (z t ), Q,^h' 2 \ from the five year Wilkinson Microwave Anisotropy Probe (WMAP) 
observations, and SDSS measurement of the baryon acoustic oscillation (BAO) scale, assuming the 
HST prior of Ho = 72 ± 8 (km/s)Mpc _1 , and without assuming spatial flatness. We find that for 
a dark energy equation of state linear in the cosmic scale factor a, the correlation of (too, too. 5) 
[too = uix{z — 0), too. 5 = vjx{z = 0.5), with wxifl) = 3too.s — 2wq + 3(wo — too. 5) a] is significantly 
smaller than that of (too, w a ) [with wx(a) = too + (1 — a)w a ]. In order to obtain model-independent 
constraints on dark energy, we parametrize the dark energy density function X(z) = px(z)/ px{0) as 
a free function with X0.5, X1.0, and X1.5 [values of X(z) at z = 0.5, 1.0, and 1.5] as free parameters 
estimated from data. If one assumes a linear dark energy equation of state, current observational 
data are consistent with a cosmological constant at 68% C.L. If one assumes X(z) to be a free 
function parametrized by (X0.5, A1.0, A1.5), current data deviate from a cosmological constant at 
z = 1 at 68% C.L., but are consistent with a cosmological constant at 95% C.L.. Future dark 
energy experiments will allow us to dramatically increase the FoM of constraints on (too, too. 5), and 
of (X0.5, A1.0, A1.5). This will significantly shrink the dark energy parameter space to either enable 
the discovery of dark energy evolution, or the conclusive evidence for a cosmological constant. 

PACS numbers: 98.80.Es,98.80.-k,98.80.Jk 
Keywords: Cosmology 



I. INTRODUCTION 

The understanding of dark energy, the unknown cause 
for the observed cosmic acceleration [l|, Hjj, continues 
to be one of the most important challenges in cosmol- 
ogy today. Dark energy could be an unknown energy 
component [H, 0, [B, B, "M Bl, or a modification of gen- 
eral relativity % Hof [n|, [ll El- and [H| con- 
tain reviews of many models. Much work continues 
to be done on the theoretical front, see for example, 
[3 El H El, US HH. Current observational data do 
not provide stringent constraints on dark energy, and 
allow a wide range of possibilities including dark en- 
ergy _being_ a cosmological constant (see, for exam ple , 




Future dark energy experiments that are significantly 
more ambitious than current ones are required to illumi- 
nate the nature of dark energy. In order to compare pro- 
posed future dark energy experiments in a useful manner, 
we need to choose the appropriate figure of merit (FoM) 
for dark energy constraints [50j ]. 

In this paper, we explore the optimization of FoM us- 
ing current observational data from supernovae, galaxy 
clustering, and cosmic microwave background anisotropy 
(CMB) data. We describe our method in Sec. II, present 



our results in Sec. Ill, and conclude in Sec. IV. 



II. METHOD 

A. General definition for figure of merit 

When we estimate a set of parameters, {fa} (i=l, 2, 
N), from data, the most intuitive figure of merit is the 
A^-dimensional volume enclosed by the 68% or 95% confi- 
dence level (C.L.) contours of the parameters. If the like- 
lihood surfaces for all the parameters are Gaussian, this 
./V-dimensional volume is proportional to the squre root 
of the covariance matrix of {fa}, ydet Cov(/i, fa, fa, ...). 
For N = 2, the 68% or 95% C.L. contours of fa and fa are 
ellipses, with enclosed area given by -Ky/det Cov(/i, fa) 
multiplied by 2.30 or 6.17. Parametrizing dark energy 
equation of state as ioj(o) = u>o + (1 — a)w a [5l[, the 
Dark Energy Task Force (DETF) defined FoM to be the 
inverse of the area enclosed by the 95% C.L. contour of 
(w , w a ) [52J, i.e., 



FoMdetf = 



1 



6. 17tt a (w a) c{w p ) 



(1) 



where w p — wq — w a (SwoSw a ) / (5w^) , and a(wi) — 
\J (Swf). Note that a(w a )a(w p ) = \J det Cov (w ,w a ), 
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thus the conversion to w p is not needed to calculate the 
FoM. 

For real data, the likelihood surfaces for the param- 
eters {fi} are almost always non- Gaussian at the 95% 
C.L., thus defining the FoM as enclosed area or volume 
by the 95% C.L. contours of {fi} becomes problematic. 
We propose the definition for a relative generalized FoM 
given by 



FoM r 



VdetCov(A,/ 2 ,/ 3 ,...r 



(2) 



where {fi} are the chosen set of dark energy parame- 
ters. This definition has the advantage of being easy 
to calculate for either real or simulated data. We have 
streamlined the definition to omit numerical factors since 
what matters is the relative FoM between different exper- 
iments. 

Note that while this FoM definition has an intuitive 
physical interpretation, it rewards experiments that yield 
very correlated estimates of the dark energy parameters. 
This is especially true in applying the DETF FoM, since 
(wQ,w a ) are always highly correlated. Hence the dark 
energy FoM [as defined in Eq.J2|)] is most meaningful 
when the dark energy parameters {fi} are chosen such 
that they are minimally correlated with each other. 



B. Dark energy parametrization 

We study constraints on a 2-parameter dark energy 
equation of state Wx(z) linear in a, as well the dark 
energy density function X(z) = px(z)/ px(0) as a free 
function at z < 1.5. 

The 2-parameter wx (z) is given by 



w x {a) 



w 



1 



a c — 1 J \a c — 1 

dcWo — w c + a(w c — wq 
a r — 1 



w c 



(3) 



where wq — wx{z — 0), and w c — wx(z — z c ). Eq.© 
corresponds to a dark energy density function 



X{z) — exp < 3 



1 



w c 



a c w - w c 

a c -l 
-w Q \ ; 



ln(l + z)+ 



(4) 



i c — 1 J 1 + z 
Eq.Q is related to wx(z) = wq + (1 — a)w a by setting 
w c - w a 



1 



or 



w c = w + (l-a c )w a . (5) 



If we choose a c = 1 + a 2 (wo)/u 2 (woW a ), then (wo,w c ) 
are uncorrelated. For current data, z c ~ 0.3. Choosing 
a c to make (wo,w c ) uncorrelated has the disadvantage 
that a c is different for different data sets. 



We recommend choosing z c = 0.5; it is sufficiently close 
to z c ~ 0.3 that the correlation of Wo and W0.5 = wx(z — 
0.5) is relatively small. It is straightforward to show that 
if \a 2 (wow a ) /[o~(wo)a(w a )}\ < 1, (wo,w c ) are less corre- 
lated than (wo,w a ) if 

<J 2 {wo) < 2 |(1 - a c )a 2 (w w a )\ ■ (6) 

This is always satisfied for z c = 0.5. Choosing z c = 
0.5, the correlation of (wo,wo.5) is smaller than that of 
(wo,w a ) by about a factor of 2 for the combined SNe, 
BAO, and CMB data considered in this paper. Fixing z c 
has the significant advantage of allowing the comparison 
of the same dark energy property for different data sets. 
For our results for the 2-parameter dark energy model, 
we use Eq.© with a c = 1/(1 + 0.5) = 2/3 (z c = 0.5). 
Thus 



w x {a) = 3w .5 - 2w + 3 (w - too.s) a 



(7) 



X(Z) = (1 + z )3(l-2» +3»o.5) exp 



9(w ~ wo. 5) 



1 + z 



In order to obtain model-independent constraints on 
dark energy, we parametrize the dark energy density 
function X(z) = px(z)/ px(0) as a free functions with 
^0.5) -X~1.0) and X1.5 [values of X(z) at z — 0.5, 1.0, and 
1.5] as free parameters estimated from data. At z > 1.5, 
we choose either X(z) = Xi. 5 , or X(z) = X x . 5 e '^ 1 - 51 
(with a as an additional parameter to be estimated from 
data). Our results are insensitive to the assumption 
about X(z) at z > 1.5 (other than that dark energy 
becomes insignificant at z > 1.5). As more data become 
available at z > 1.5, we can include X2.0 = X(z = 2.0), 
X2.5 = X(z = 2.5), and X3.0 = X(z = 3.0) as esti- 
mated parameters, as well as inserting more estimated 
X{z) values at z < 1.5. Early dark energy (significant 
at high z) is not required by current data, and leads to 
contradiction with observed cosmic structure formation 
[69l |. unless a cutoff is imposed. 

The constraints on X0.5, X\,o 1 and X1.5 are insensitive 
to the interpolation used in deriving X(z) elsewhere. The 
simplest smooth interpolation is given by a polynomial: 



X(z) 



--(3Z-l)(3Z-2)(2-l) 
+ ^X Q . 5 z(3z-2) (z-1) 



A!. z(3z-l)(z-l) 



+-Xi. 5 2(3z-l)(32-2) ! 



(8) 



where z = z/1.5. 



C. Data analysis technique 

The comoving distance from the observer to redshift z 
is given by 

r{z) = cH^ 1 [Qkl-^sirmUQk^ 2 T(z)}, (9) 
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T(z) 



dz' 



E(z) = H(z)/H 



where 0/. = —k/Hg with k denoting the curvature con- 
stant, and sinn(x) = sin(x), x, sinh(ir) for f2/- < 0, 
fit = 0, and Qk > respectively, and 

e 2 {z) = n^i+zf+n^dii+zy+rikii+zf+nxXiz) 

(10) 

with Qx = 1 — — ^rad — ^fe, and the dark energy 
density function X(z) = px(z)/px(0)- 

CMB data give us the comoving distance to the 
photon-decoupling surface r(z„), and the comoving sound 
horizon at photo-decoupling epoch [53|, HH 



t» 



cm 



dt 



dz — r-r- 

.. E{z) 
da 



3(1 + R b a) a i E 2 {z) 



(11) 



where p^ ata are the maximum likclyhood values given in 
Table 10 of [E^. 

SN la data give the luminosity distance as a function 
of redshift, d L (z) = (1 + z)r(z). We use 182 SNe la 
from the HST/GOODS program [59| and the first year 
SNLS 60], together with nearby SN la data, as compiled 
by [59^ We do not include the ESSENCE data [HI, 
as these are not yet derived using the same method as 
thosed used in [59| . Combining SN la data derived using 
different analysis techniques leads to systematic effects 
in the estimated SN distance moduli [6l|, [62j • Appendix 
A of [HI describes in detail how we use SN la data (flux- 
averaged to reduce lensing-like systematic effects [62|, |H, 
l&ij and marginalized over Hq) in this paper. 

We also use the SDSS baryon acoustic oscillation 
(BAO) scale measurement by adding the following term 
to the x 2 °f a model: 



Xbao 



(A - Abao) 



(17) 



where a is the cosmic scale factor, a* = 1/(1 + 2*), 
and a A E 2 (z) — fi m (a + a cq ) + Q k a 2 + VL x X(z)a A , 



where A is defined as 



with 



= fj ra d/fim = 1/(1 + Zeq), and Z c 



2.5 x 10 4 ^ m fe 2 (r CM B /2.7K)- 4 . The sound speed is 

c s = l/y3(l + Rb a), with R b a = 3 i o f) /(4p 7 ), R b = 
31500n b h 2 (T C MB/2.7K)- 4 . COBE four year data give 
Tomb = 2.728 ± 0.004 K (95% C.L.) [Hi]. We take 
T C mb = 2.725 following pftt, since we will use the CMB 
bounds derived by [E^]- The angular scale of the sound 
horizon at recombination is defined as l a — 7rr(z*)/r s (z*) 

Wang & Mukherjee 2007 [5J] showed that the CMB 
shift parameters 

R = \Jn m HQr(z*), l a = TTr(z„)/r s (z^) : (12) 

together with fl b h 2 , provide an efficient summary of 
CMB data as far as dark energy constraints go. We 
use the covariance matrix of [R(z*), l a (z*), £l b h 2 ] from the 
five year WMAP data (Table 11 of ^3), with z^given 
by fitting formulae from Hu & Sugiyama (1996) [581 ]: 

z* = 1048 [1 + 0.00124(fM 2 )-°- 738 ] [l + gi {Q m h 2 ) 92 ] , 

(13) 

where 



-0.238 



,9i 



92 = 



0.0783 {fl b h 2 ) 
1 + 39.5 (n b h 2 ) - 763 

0.560 
l + 21.1(fi 6 /i 2 ) 1 - 81 



(14) 
(15) 



CMB data are included in our analysis by adding the 
following term to the x 2 of a given model with pi — 
R(z*), p 2 = l a {z*), and p 3 = fl b h 2 : 



XCMB 



Api [C ov 1 (pi ,pj)] Ap 3 , Ap l =pi-p 



data 



(16) 



r 2 (z B Ao) 



czbao 
H(zbao) 



1/3 



(O r 



Hi 



1/2 



CZBAO 



(18) 



and Abao = 0.469 (n s /0.98)-° 35 , a A = 0.017, and 
zbao = 0.35 (independent of a dark energy model) [65] ]. 
We take the scalar spectral index n,<? = 0.96 as measured 
by WMAP five year observations [571 ] . 

For Gaussian distributed measurements, the likelihood 



function L 



oc e 



-x 2 / 2 , with 



2 _ 2 
X — XCMB 



X 2 SNc 



Xbaoi 



(19) 



where Xcmb * s given in Eq. p^)) . XsNe ^ s given in Ap- 
pendix A of [56(, and Xbao i s given in Eq. pT|l . 

The current BBN constraints 66], S = 0.942 ± 0.030 
{N v = 2.30lo;|!) rule out the standard model of particle 
physics (S — 1, N v — 3) at la [66[. Given the uncertain- 
ties involved in deriving the BBN constraints, we relax 
the standard deviation of S by a factor of two, so that 
the standard model of particle physics is allowed at la. 
We find that the resultant BBN constraints do not have 
measurable effect on our dark energy constraints. 

For all the dark energy constraints from combining the 
different data sets presented in this paper, we marginalize 
the SN la data over Hq in flux-averaging statistics [56| . 
and impose a prior of Hq = 72 ± 8 (km/s)Mpc~ 1 from 
the HST Cepheid variable star observations [671 ] . 

We run a Monte Carlo Markov Chain (MCMC) based 
on the MCMC engine of (68| to obtain O(10 6 ) samples 
for each set of results presented in this paper. The pa- 
rameters used are (f2fc, f2 TO , h, fl b h 2 , Pbe)- The dark 
energy parameter set is described in Sec.IIB. We as- 
sumed flat priors for all the parameters, and allowed 
ranges of the parameters wide enough such that further 
increasing the allowed ranges has no impact on the re- 
sults. The chains typically have worst e-values (the vari- 
ance(mean)/mean(variance) of 1/2 chains) much smaller 
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than 0.005, indicating convergence. The chains are sub- 
sequently appropriately thinned to ensure independent 
samples. 



III. RESULTS 

FigQ] shows the 68% and 95% C.L. contours of 
(wojWo.s) (upper panel) and (wo,w a ) (lower panel) from 
WMAP 5 year measurement of [R(z*), l a (z*), Qbh 2 ], 
and 182 SNe la (from the HST/GOODS program, the 
first year Supernova Legacy Survey, and nearby SN la 
surveys), with and without the SDSS measurement of 
the baryon acoustic oscillation (BAO) scale. We have 
assumed the HST prior of H = 72 ± 8 (km/s)Mpc" 1 , 
and allowed Qk to vary. Tabic 1 shows the mean, rms 
variance, and correlation coefficients of (100,100.5) an d 
(wo,w a ), as well as the relative dark energy FoM r de- 
fined in Eq.0. Note that Pearson's correlation coeffi- 
cient p xy = a 2 (xy)/[a(x)a(y)}. Adding the SDSS BAO 
scale measurement improves the FoM r by a factor of 21.5 
for (W07W0.5), and by a factor of 27.0 for (wo,w a ). Since 
(wo,wo.5) are significantly less correlated than (wo,w a ), 
the improvement factor in FoM r of (wojifo.s) is a more 
reliable indicator of the impact of adding the SDSS BAO 
scale measurement. 

Fig. 2 shows the one dimensional marginalized prob- 
ability distributions (pdf) of (fi m , h, Clk, fl b h 2 , Wq, 
w a ), for 182 SNe la, the SDSS BAO scale measurement, 
and the WMAP 5 year data in the form of measured 
(1) [R(zJ, l a M, tt b h 2 } (solid lines), (2) [i?(z»), 
il b h 2 } with z* fixed at 1090.4 (dotted), and (3) [R{zj, 
l a {z*), z*] (dashed). For reference, the dot-dashed line 
shows the pdfs for 182 SNe la, the SDSS measurement of 
the BAO scale, and the WMAP 3 year data in the form of 
measured \R(zqmb), l a (z C MB), SI&/1 2 ] with z C mb fixed 
at 1089 from [56[. 

We find that in spite of the different pdfs for fl b h 2 , 
using the [R(z*), l a (z*), Q b h 2 } and [R(z*), l a (z*), 2*] 
measurements give about the same constraints on (f2 m , 
h, Qk, w , w a ). Using the [R(z*), l a (z*), Q b h 2 ] mea- 
surement with z* fixed at 1090.4 gives slightly tighter 
constraints on (wq, w a ). In combination with the super- 
nova and BAO data, the WMAP 5 year data improve 
constraints on (wq, w a ) slightly compared to the WMAP 
3 year data, while tightening the constraints on Qk and 
h. 

Fig. 3 shows the constraints on the dark energy density 
function X(z) — px(z)/ px(0) parametrized by Eq.([5]), 
with X(z) at z > 1.5 given by either X(z) = A1.5 or 
X(z) — Xi,5 exp[a(z — 1.5)]. Note that the assump- 
tion about dark energy at z > 1.5 has only a weak ef- 
fect on the dark energy constraints at z < 1.5. Note 
that taking X(z > 1.5) = X1.5 gives slightly less strin- 
gent constraints on dark energy at z < 1.5. This is 
because parametrizing dark energy at z > 1.5 with an 
extra parameter requires choosing the early dark energy 
parametrization such that it is not degenerate with cos- 
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FIG. 1: The 68% and 95% C.L. contours of (u>0i u '0.5) (upper 
panel) and (wo, w a ) (lower panel) from WMAP 5 year measurement 
of [R(z m ), L{z*), n b h 2 }, and 182 SNe la (from the HST/GOODS 
program, the first year Supernova Legacy Survey, and nearby SN la 
surveys), with and without the SDSS measurement of the baryon 
acoustic oscillation (BAO) scale. 



mic curvature; this is why f2fe is not well constrained if 
we choose X(z > 1.5) = X1.5 (l + z) a , but £lk is well con- 
strained if we choose X(z > 1.5) = A1.5 exp[a(z — 1.5)] 
[56| . The latter helps break the degeneracy oi with 
X(z), thus leading to much tighter constraints on f2/- 
and slightly tighter constraints on X(z) at z < 1.5 (see 
Fig. 3). This suggests that the more conservative ap- 
proach in constraining dark energy is to assume that 
X(z>1.5)=X 1 . 5 . 

Table 2 shows the mean, rms variance, and correlation 
coefficients of (X 5 , X\ .0, X 15 ), as well as the relative 
dark energy FoM r defined in Eq.([2]). 
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TABLE I: Constraints on (two, W0.5) and (wo,w a ) 



Data 


fi(w ) 


a(w ) 






1 P™0™0.5 


FoM r 


WMAP5+SNe 


-1.075 


0.598 


-1.939 


1.572 


-0.401 


1.163 


WMAP5+SNe+BAO 


-0.937 


0.226 


-0.953 


0.206 


-0.512 


25.013 


Data 


fi(w ) 


a(w ) 




0-(w a ) 


pWQW a 


FoM r 


WMAP5+SNe 


-1.073 


0.647 


-2.960 


6.759 


-0.670 


0.308 


WMAP5+SNe+BAO 


-0.938 


0.226 


-0.045 


1.126 


-0.882 


8.326 



TABLE II: Constraints on X(z) parametrized by (A0.5, X1.0, X1.5) 



X(z > 1.5) n(Xo.s) cr(X . 5 ) fx(Xi. ) cr(Xi.o) /u(^i.s) er(Jfi. 5 ) px 5 x L0 px 5 x L5 px 1() x L5 FoM r 
Xi.s 1.059 0.213 2.556 1.215 7.503 8.037 -0.389 -0.666 0.906 2.0771 

Xi. 5 e a(z " u) 1.091 0.195 2.436 1.121 6.533 7.351 -0.303 -0.609 0.895 2.402 
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FIG. 2: One dimensional marginalized pdfs of (f2 m , h, Q^, Q b h 2 , 
wo, w a ) from 182 SNe la, the SDSS BAO scale measurement, and 
the WMAP 5 year data in the form of measured (1) [R(z„), i a (z»), 
Q b h 2 ] (solid lines), (2) [_R(z»), l a (z*), Q b h 2 ] with 2* fixed at 1090.4 
(dotted), and (3) [i?(z*), £ a (z„), z»] (dashed). The dot-dashed 
line shows the pdfs for 182 SNe la, the SDSS measurement of the 
BAO scale, and the WMAP 3 year data in the form of measured 
[R(zcmb), Ici(zcmb), 0. b h 2 ] with zcmb fixed at 1089 from pjjj. 



IV. SUMMARY AND DISCUSSION 

In order to compare current and future dark energy 
experiments on the same footing, we have introduced a 
simple and straightforward definition for the Figure-of- 
Merit (FoM) of constraints on any set of dark energy 



FIG. 3: One dimensional marginalized pdfs of dark energy and cos- 
mological parameters from WMAP 5 year measurement of [i?(z*), 
l a (z*), fl b h 2 ], 182 SNe la, and the SDSS BAO scale measure- 
ment. Solid and dashed lines indicate X(z) at z > 1.5 given by 
X(z) = Xi.s and X(z) = X1.5 exp[ct(z — 1.5)] respectively. 



parameters, Eq.([2]), that is easily applicable to both real 
and simulated data. 

We recommend the adoption of two dark energy 
parametrizations in comparing different experiments: (1) 
A dark energy equation of state wx(z) linear in a, with 
its values at z = and z = 0.5, (wo, W0.5), as param- 
eters estimated from data [see Eq. (|7j)] . We find that 
(w 0l w 5) are significantly less correlated than (wo, w a ) 
[see Table 1 and FigJT], hence the factor of improve- 
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ment in the FoM r [as defined in Eq.j2j] for (wq, iuo.s) 
is a more reliable indicator of the improvement in dark 
energy constraints than the factor of improvement of 
FoM r for (wq, w a ). (2) The dark energy density func- 
tion X(z) = px(z)/ px(0) parametrized by its values at 
z = 0.5, 1.0, and 1.5, (X0.5, X1.0, X1.5), for z < 1.5 
[see Eq.©], and X(z > 1.5) = Xt, 5 . We find that this 
flat cutoff in X(z) gives more conservative constraints on 
X(z) than parametrizing early dark energy with an extra 
parameter such that cosmic curvature is constrained (see 
Fig®. 

We have demonstrated the use of the FoM r [see Eq.©] 
for these two dark energy parametrizations [sec Eq.(|Tj) 
and Eq.flHJ)] using WMAP 5 year measurement of [R(z*), 
l a (z*), Slbh 2 ], 182 SNe la (from the HST/GOODS pro- 
gram, the first year Supernova Legacy Survey, and nearby 
SN la surveys), and the SDSS measurement of the baryon 
acoustic oscillation (BAO) scale [see FigsJTJS] . Dark en- 
ergy is consistent with a cosmological constant at 68% 
C.L. if one assumes the two-parameter dark energy equa- 
tion of state model, wx(a) = 3^0.5 — 2wo + 3(wq — w 5 ) a. 
If one assumes dark energy density to be a free func- 
tion parametrized by its values at z = 0.5, 1.0, and 1.5, 
then dark energy deviates from a cosmological constant 
at z = 1.0 at 68% C.L., but is consistent with a cosmolog- 
ical constant at 95% C.L. (see FigJ3]). This illustrates the 
importance of using the model-independent parametriza- 
tion in probing dark energy. Measuring X(z) as a free 
function from data allows us to detect epochs of vari- 
ation in dark energy density. It also allows us to con- 



strain a broader class of dark energy models than rep- 
resented by wx{z); for example, dark energy models in 
which X{z) becomes negative in the past or future, which 
are excluded by fiat if one only measures Wx(z) since 
X{z) = exp{f* dz' 3[1 + w x (z'j}/(l + z')} [H[. The two 
parameter dark energy equation of state model (linear in 
a) implies strong assumptions about dark energy, and is 
not sensitive to a transient variation in dark energy; thus 
it is most useful in comparing forecasts for future dark 
energy experiments under the simplest assumptions. 

Future dark energy experiments from both ground and 
space [H Izljll IZa IZa Il3, together with CMB data 
from Planck [78j, will dramatically improve our ability 
to probe dark energy, and eventually shed light on the 
nature of dark energy. Using both a linear dark en- 
ergy equation of state [parameterized by (wo, W0.5)} and 
dark energy density function X(z) as a free function 
[parametrized by (X0.5, -X1.0, -X1.5)] provides a simple 
and balanced approach to exploring dark energy. Pro- 
posed future dark energy experiments should be evalu- 
ated by comparing their FoM r for both (wq, ifo.s) and 
(X0.5, Xt,o, Xl.s) to that of current data. 
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